
* counts from counterfactual *
sum pcount if ltd_iss_all_other > 9500 & ltd_iss_all_other <= 10000 
scalar cfbcount = r(sum)
sum pcount if ltd_iss_all_other > 10000 & ltd_iss_all_other <= 14500 
scalar cfacount  = r(sum)

* actual counts  * 
sum count if ltd_iss_all_other > 9500 & ltd_iss_all_other <= 10000 
scalar actbcount = r(sum)
sum count if ltd_iss_all_other > 10000 & ltd_iss_all_other <= 14500 
scalar actacount = r(sum)

* size of excess mass
return scalar totmanip_ex = (actbcount - cfbcount) / totalcount 
return scalar inranmanip_ex = (actbcount - cfbcount) / cfacount 

* size of missing mass
return scalar totmanip = (cfacount - actacount) / totalcount 
return scalar inranmanip = (cfacount - actacount) / cfacount